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Abstract 

The difference between the ratio of multiplicities and the ratio 
of their derivatives on energy for gluon and quark jets is calculated 
up to next-to-next-to leading order of perturbative QCD. Its non¬ 
zero value is uniquely defined by the running property of the QCD 
coupling constant. It is shown that this difference is rather small 
compared to values which can be obtained from experimental data. 

This disagreement can be ascribed either to strong non-perturbative 
terms or to experimental problems with a scale choice, jets separation 
and inadequate assignement of soft particles to jets. 

The ratio of multiplicities in gluon and quark jets has been of much debate 
during last years. The lowest order prediction of the perturbative QCD dras¬ 
tically overvalues experimental results for this ratio. Therefore, the higher- 
order corrections up to the next-to-next-to leading terms of the perturbative 
QCD have been calculated §0. and the ratio has been written as 

= r (v) = £r(! - D7o (y) - wjfty)) + O(7o)> 4) 


where the scale of the process is defined by y — lnQ/Q 0 , Q is a virtuality of 
a jet, Q 0 =cons t; 
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Here, rif is the number of active flavours, N c =3 the number of colours. Let 
us note that only the first term of r 2 was obtained in |l| by Feynman graphs 
technique. Other terms come from higher derivatives of the generating func¬ 
tions || when the solution of the QCD equations for generating functions 
is found out by Taylor series expansion. They take into account the energy 
conservation in three-parton vertices. 

Asymptotically, the ratio r(y) (jlj) tends to a constant Ca/Cf= 9/4 which 
shows that the gluon jet bremsstrahlung is stronger than for quark jets |||. 
However, the correction terms are still noticeable at presently accessible en¬ 
ergies. 

The energy dependence of average multiplicities of gluon and quark jets 
used to be expressed in terms of the anomalous dimension 7 (y) as 


(n G (y)) = Ke$ 7(2/W , (n F (y)) = (n G (y))/r{y). ( 6 ) 


Then it is easy to get the ratio of their derivatives as 
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or for its difference with the ratio of multiplicities one gets 
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The anomalous dimension in terms of the running coupling constant looks 
like 

7 = 7o(l-ai7o) + 0(7o). (9) 
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The higher order terms have been omitted in our treatment, 
account that 


7o = -hi7o + O(7o), 
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one estimates that the second term in the denominator of ( 0 ) can be neglected 
and gets finally 

+ ( 12 ) 

The formulas (§) and (0 are the main result of the paper. 

It is important to stress that the difference D(y) between the ratios of 
derivatives and multiplicities directly demonstrates the running property of 
the QCD coupling constant. It is identically equal to zero for fixed cou¬ 
pling because r' — 0 in that case. For the running coupling, this difference 
is always positive and proportional to the value of the coupling constant. 
Consequently, it tends to zero in asymptotics. 

The corrections due to the anomalous dimension (0) and due to the multi¬ 
plicity ratio ( 0 ) have been specially left as separate factors in brackets in (|T^) 
to show their relative importance. One easily notices that the last correction 
in D(y ) is stronger than in the ratio (0) itself because r' in the numerators of 
( 0 ) and ( 0 ) acquires the factors n when differentiating the subsequent terms 
of 7q in (0). 

Asymptotically the terms in brackets tend to 1 since y 0 —»► 0, but at 
present energies (70 ~ 0.45 — 0.5 at Z°) the corrections are rather large. 
In the energy range near Z°, the factor in front of the brackets in (|l2|) is 
about 0.05-0.06. The anomalous dimension correction in the first bracket 
enlarges the value of D by about 15%. More important is the expression in 
the second bracket. With values of 77 and r 2 given by (3!), (0) and inserted 
in (|T2|), one estimates it as about 3.2, i.e. the linear in y 0 term (NNLO 
correction) contributes more than 1 (MLLA term), and therefore next terms 
should be evaluated. Anyway taking this expression for granted, one would 
estimate D as 

D « 0.16 — 0.20. (13) 

Thus the predicted value of D is comparatively small. 

Let us note that the first term from (0) only as calculated first in [|TJ 
contributes to the second bracket the value about 0.25, and the corresponding 
contribution to D is equal approximately 0.06 - 0.08. 

We have considered the one-scale problem without imposing any limita¬ 
tions on the transverse momentum and considering the jet evolution with 
energy in the total phase space. In experiment, the more severe restrictions 
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are often introduced due to specifics of installations or due to some special 
criteria. In particular, one can consider jet evolution as a function of some 
internal scale (e.g., similar to the jet transverse momentum) at a given energy 
(usually chosen at Z° peak because of larger statistics). Such a procedure has 
been exploited in Refs |j. It was claimed that the proportionality of the 
two scales favours comparison of theoretical predictions with experimental 
data on the evolution of above ratios. The scale chosen in || was 
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where Ej et is the energy of the jet and 9 the angle with respect to the closest 
jet. In |Q], the similar scale of the geometric mean of such scales of the gluon 
jet with respect to both quark jets in three-jet events has been used. From 
the data presented in [Q one can estimate that the difference D(k) ranges 
approximately from 0.95 at k= 7 GeV to 0.7 at k= 28 GeV. These values are 
much larger than predicted above. 

However, both theoretical and experimental approaches to the problem 
should be further analysed. The large correction in Eq. (|l^) implies that 
next terms of the perturbative expansion can be essential in the ratio of 
derivatives. They are enlarged in (|12|) as discussed above and support earlier 
conclusions (see, e.g., Ref. |J) that subleading effects are quantitatively 
important at experimentally accessible energies. Our experience tells us also 
that power corrections due to conservation laws related to recoil effects can 
become essential as well |Qj. 

Experimental procedures raise even more questions. There exists the scale 
dependence in the final result if one compares the data of || with that of 
[j5j]. The problem of the scale choice has recently been discussed in Ref. [§]. 
Besides, it has been shown || that different algorithms of the jet selection 
and ascribing soft particles to a jet give rise to somewhat different conclusions 
concerning jet multiplicities and their slopes. Further study is necessary. 

Another problem raised by the procedure used in Ref. || is the role 
of non-perturbative effects in average multiplicities of gluon and quark jets. 
It has been claimed || that the hadronization of quark jets leads to some 
constant (an assumed ansatz!) excess over gluon jets equal to about two 
or three particles. It compensates somewhat the perturbative QCD excess 
of the gluon jet multiplicities. This constant term in energy dependence 
of multiplicities has been ascribed to a possible non-perturbative contribu- 
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tion. Unfortunately, this statement can hardly be confronted to any reliable 
theoretical treatment. These terms do not contribute to derivatives of the 
multiplicities but drastically diminish the multiplicity ratio thus increasing 
their difference D(y). 

At the same time, the theoretical estimate ( [dip seems quite reasonable 
qualitatively if one accepts the values of the ratio of derivatives as about 2.1 
given in || and the value of the ratio of multiplicities about 1.84 predicted 
theoretically || with the experimental result ~ 1.6 claimed in 0]. 

In conclusion, we have evaluated the difference between the ratio of slopes 
of average multiplicities and the ratio of multiplicities in gluon and quark jets 
in the next-to-next-to leading approximation of the perturbative QCD and 
confronted it to some experimental data. 
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